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Abstract 

Activating transcription factor 3 {Atf3) is rapidly and transiently upregulated in numerous systems, and is associated with 
various disease states. Atf3 is required for negative feedback regulation of other genes, but is itself subject to negative 
feedbacl< regulation possibly by autorepression. In cardiomyocytes, Atf3 and Egrl mRNAs are upregulated via ERKl/2 
signalling and Atf3 suppresses Egrl expression. We previously developed a mathematical model for the Atf3-Egr1 system. 
Here, we adjusted and extended the model to explore mechanisms of Atf3 feedback regulation. Introduction of an 
autorepressive loop for Atf3 tuned down its expression and inhibition of Egrl was lost, demonstrating that negative 
feedback regulation of Atf3 by Atf3 itself is implausible in this context. Experimentally, signals downstream from ERKl/2 
suppress Atf3 expression. IVlathematical modelling indicated that this cannot occur by phosphorylation of pre-existing 
inhibitory transcriptional regulators because the time delay is too short. De novo synthesis of an inhibitory transcription 
factor (ITF) with a high affinity for the Atf3 promoter could suppress Atf3 expression, but (as with the Atf3 autorepression 
loop) inhibition of Egrl was lost. Developing the model to include newly-synthesised miRNAs very efficiently terminated 
Atf3 protein expression and, with a 4-fold increase in the rate of degradation of mRNA from the mRNA/mlRNA complex, 
profiles for Atf3 mRNA, Atf3 protein and Egrl mRNA approximated to the experimental data. Combining the ITF model with 
that of the miRNA did not improve the profiles suggesting that miRNAs are likely to play a dominant role in switching off 
Atf3 expression post-induction. 
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Introduction 

The concept of immediate early genes (lEGs) was initially 
established in relation to viral infection of bacterial or eukaryotic 
cells (e.g. herpesvirus [1]). Here, lEGs are defined as the genes 
encoding the first phase of mRNAs expressed from the viral 
genome, relying entirely on pre-existing host proteins. Early and 
delayed genes are expressed later and require production of new 
proteins. In mammaUan cells, lEGs became defined as those 
which are regulated by pre-existing transcription factors (TFs) such 
that increases in lEG mRNA expression are not suppressed by 
protein synthesis inhibitors (e.g. cycloheximide) [2]. Many lEGs 
encode transcriptional regulators that are required to modulate 
downstream gene expression (i.e. second phase genes) and some 
are inhibitory factors required to terminate transcription. With the 
discovery of miRNAs that influence mRNA expression and 
translation [3,4], some adjustment of the lEG concept is necessary 
according to whether mRNA regulation is transcriptional or post- 
transcriptional. Synthesis, production and degradation of miRNAs 
may be insensitive to inhibitors of protein synthesis as would be 
any (post-transcriptional) effects on mRNA expression. 



Cardiomyocytes, the contractile cells of the heart, withdraw 
from the cell cycle perinatally. Neonatal rat cardiomyocytes are 
therefore highly synchronized and, because they do not divide, 
form a good model for the study of lEG regulation. In these cells, 
endothelin-1 (ET-1), a Gq protein-coupled receptor (GqPCR) 
agonist, elicits maximal activation of the entire pool of extracel- 
lular signal-regulated kinases 1/2 (ERKl/2, the prototypic 
mitogen-activated protein kinases) within 3-5 min [5,6]. The 
earliest lEGs are substantially and maximally upregulated within 
15-30 min, with second phase RNAs being detected within 1 h 
[7]. ERKl/2 play a major role in the response and promote 
upregulation of ~70% of ET- 1 -responsive transcripts [7,8]. 
Recent data indicate that ERKl/2 act together with their 
downstream substrates p90 ribosomal S6 kinases (RSKs) to 
regulate RNA expression, with RSKs being required to increase 
expression of ~50% of the RNAs upregulated by ET-1 [9]. As in 
other systems, upregulation of cardiomyocyte lEG mRNAs is 
acute and transient. For example, activating transcription factor 3 
(Atf3) and early growth response 1 (Egrl) mRNAs (each regulated 
via ERKl/2) are both upregulated maximally by ET-1 within 
~30 min [10]. However, expression of Egrl mRNA returns to 
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Author Summary 

Activating transcription factor 3 (Atf3) is an important 
regulatory transcription factor which is associated with 
inflammation, restraint of the immune response and 
cancer. In this work, we develop a series of mathematical 
models to understand how Atf3 may be regulated. 
Informed with data from the literature and our own 
experiments, we show that self-regulation of Atf3 does not 
allow for variation between experimentally observed Atf3 
mRNA and Atf3 protein expression profiles. A fast-acting 
signal via phosphorylated RSK is also shown to be 
implausible for similar reasons. Extending our mathemat- 
ical model further, we postulate for the first time, that the 
observed dynamical variation in Atf3 mRNA and protein 
can be described by microRNAs downstream of RSKs. The 
further inclusion of an inhibitory transcription factor for 
Atf3 expression has little effect on these findings. 

basal by ~2 h and expression of Atf3 mRNA declines close to 
basal levels over ~4 h. This raises the question of negative 
feedback regulation. 

Atf3 is emerging as an extremely important feedback regulator 
of transcription. Particular emphasis is placed on its role in 
inflammation and Atf3 is essential for restraining the immune 
response, but AtB is also associated with cancer and cardiac 
dysfunction [11-14]. However, although it is stress-regulated and 
stress-regulating, Atf3 is upregulated in many systems by growth 
stimuli including peptide growth factors and GqPCR agonists 
[7,15-17]. Atfi forms homo- or heterodimers that bind to ATF/ 
CRE sites in gene promoters to regulate transcription [18]. It is 
largely regulated at the level of expression, being present at very 
low levels in quiescent cells and induced as an lEG by a range of 
extracellular stimuli and cellular stresses. ERKl/2 are particularly 
implicated in promoting Atf3 mRNA expression and several TFs 
may be involved [19-21]. AtB is generally viewed as a 
transcriptional repressor, particularly when acting as homodimers, 
and it may repress transcription from its own promoter to limit 
expression [22]. In cardiomyocytes, Atf3 operates in a negative 
feedback system with Egrl and, by binding to the Egrl promoter, 
AtB inhibits Egrl transcription [10]. Mathematical modelling of 
the system demonstrated that this, in itself, could suffice for the 
transient nature of the Egrl response. The pivotal role that AtB 
plays in biological systems renders it essential to understand the 
mechanisms that regulate AtB expression. Here, we have 
developed our original mathematical model to determine whether 
AtB serves as a negative feedback regulator of its own 
transcription (as suggested [22]) and to explore other mechanisms 
that may switch off AtB mRNA expression. We demonstrate that 
self-regulation by AtB on its own promoter (or, indeed, self- 
regulation by any transcription factor) is implausible and miRNAs 
are likely to play a dominant role in switching off AtB expression 
post-induction. 

Results 

AtB mRNA and protein are rapidly and transiendy upregulated 
in cardiomyocytes exposed to ET-1 (Figure lA) and, since loss of 
AtB produces sustained expression of Egrl mRNA (Figure IB), 
AtB negatively regulates Egrl expression [10]. The initial 
deterministic ordinary dilferential equation model generated for 
this AtB-Egrl negative feedback system (detailed in [10]) made the 
following assumptions: AtB and Egrl are co-regulated through 
ERKl/2; ERKl/2 increase transactivating activities of TFs 



pre-bound to each promoter (i.e. the TF/promoter constitutes a 
single entity regulated by ERKl/2); since ERKl/2 activity 
remains high for at least 30 min [23], we did not switch off the 
protein kinase signal. This model predicted that AtB alone could 
suffice to to switch off Egr 1 transcription. 

To address questions relating to AtB transcription we developed 
the original model further (Figure IC). Full details of reaction 
equations, mathematical modelling and additional parameterisa- 
tion are provided in Text SI. We allowed for the rate of 
dephosphorylation (i.e. inactivation) of ERKl/2 following stimu- 
lation with ET-1, considering that phosphorylation of ERKl/2 is 
maximal at ~3 min (Figure ID, [9]) and declines by ~90% by 1 h 
with a return to baseline by 2 h [23]. The modelling profile for 
ERKl/2 phosphorylation is shown (Figure IE). We allowed for 
the rate of phosphorylation and dephosphorylation of TFs bound 
to each of the AtB and Egrl promoters. In other cells, 
phosphorylation of CREB and/ or c-Jun TFs is implicated in the 
upregulation of AtB transcription [20,24,25]. In cardiomyocytes 
exposed to ET-1, accumulation of phosphorylated (i.e. activated) 
CREB and c-Jun TFs each requires ERKl/2 signalling [26,27], 
making them prime exemplary candidates for the model. 
Experimentally, maximal phosphorylation of CREB is detected 
at ~ 10 min and declines to ~25% maximal by ~ 60 min 
(Figwe IF). The time course for phosphorylated c-Jun is delayed 
(maximal at ~ 45-60 min, declining to ~50% maximal at 2 h 
(Figure IG), because of the additional increase in expression of c- 
Juii protein [26]. ET-1 also increases expression of Myc mRNA 
(maximal at 1 h) [9] that promotes expression of AtB over a longer 
time frame in cells treated with serum [19], so the signal could be 
perpetuated over a much longer time period. Here, for the AtB 
promoter, the model was initially developed using the profile for 
phosphorylated CREB that wiU give the most rapid rate at which 
AtB transcription is terminated (Figure IH). The Egrl promoter is 
most likely regulated by ternary complexes of Elkl/SRF that can 
be stimulated either by phosphorylation of Elkl [28] or an 
increase in SRF binding [29]. Elkl is rapidly and transientiy 
phosphorylated in cardiomyocytes exposed to phenylephrine with 
a similar profile to CREB following exposure to ET-1 [30]. SRF is 
upregulated by ET- 1 with maximal expression of SRF mRNA at 
~1 h [7] and this is likely to perpetuate the response given that 
our data indicate that, in the absence of AtB expression, Egrl 
expression is sustained over at least 2 h (Figure IB). Thus, for the 
TF regulating Egrl transcription, we used a similar rate of TF 
phosphorylation as for CREB with delayed dephosphorylation. 
Binding of the TF to the Egr 1 promoter was competitive with AtB 
and was modelled in the absence or presence of AtB (Figure II). 
Egrl mRNA expression remained sensitive to negative feedback 
regulation by AtB (Figure IJ). With these adjustments, AtB 
mRNA and protein were predicted to decline (Figure IK), but the 
profiles were delayed relative to the experimental data (Figure lA). 
Thus, merely switching off the signal to promote AtB transcrip- 
tion/translation alone cannot account for its subsequent down- 
regulation, indicating that a negative feedback system also operates 
for AtB. 

Model extension 1: Atf3 protein inhibits Atf3 
transcription 

Previous studies used a protein overexpression approach with 
reporter assays for the AtB promoter to provide evidence that AtB 
operates in an autorepression loop and, by binding to an element 
immediately downstream of the TATA box, it inhibits its own 
transcription [22]. Since it is difficult to test this experimentally 
with an endogenous system, we extended the mathematical model 
developed for the AtB-Egrl feedback system to determine 
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Figure 1. Experimental data for Atf3 and Egr1 regulation in cardiomyocytes exposed to ET-1 and modifications to the original 
mathematical model for the Atf3/Egr1 negative feedback system. A, Experimental data for Atf3 mRNA (solid circles, solid line) and Atf3 
protein (open circle, dashed line) levels in cardiomyocytes exposed to 100 nM ET-1 for the times shown. Results are expressed relative to maximum 
values as means ± SEM for at least 3 independent experiments (adapted from data in [10]). B, Experimental data for Egri mRNA levels in 
cardiomyocytes exposed to 1 00 nM ET-1 for the times shown in the presence of endogenous Atf3 (open circles, dashed line) or with inhibition of Atf3 
(solid circles, solid line). Results are expressed relative to maximum values as means ± SEM for at least 3 independent experiments (adapted from the 
data in [1 0]). C, The original mathematical model [1 0] was revised (dotted lines) to account for dephosphorylation of ERK1 II and to include reversible 
binding of phosphorylated transcription factors to the Egri and Atf3 promoters. Lines with arrows indicate protein activation, DNA transcription or 
protein translation. Blocked lines indicate inhibition. Protein and mRNA degradation is indicated by 0. D, Experimental data for ERK1/2 
phosphorylation in cardiomyocytes exposed to ET-1 for the times shown. Data are means for n = 3 independent experiments and are adapted from 
the data in [9]. E, Predicted profile from the model for ERKl/2 phosphorylation. F, G, Experimental data for phosphorylation of CREB (F) and c-Jun (G) 
in cardiomyocytes exposed to ET-1 for the times shown. Data are means for at least 3 independent experiments and are adapted from data in [26,27]. 
H, The predicted profile from the model for phosphorylation of the TF bound to the Atf3 promoter was based on experimental data for CREB (F). I, 
Predicted profile from the model for phosphorylation of the TF bound to the Egri promoter was based on experimental data for CREB for the initial 
increase, but the rate of decreased phosphorylation was adjusted such that, in the absence of Atf3, the profile approximated to that of the 
experimental data (B). Profiles are shown in the presence (open circle, dashed line) or absence (solid circles, solid line) of Atf3 feedback inhibition. J, 
Predicted profiles from the model for Egri mRNA expression in the presence (open circle, dashed line) or absence (solid circles, solid line) of Atf3 
feedback inhibition. K, Predicted profiles from the model for Atf3 mRNA (solid line) or protein (dashed line) expression. 
dor:l 0.1 371/journal.pcbi.l 003597.g001 



whether a direct Atf3 autorepression loop is feasible (Figure 2A). 
The development of this and sub.sequent mathematical models is 
detailed in Text SI. The model extension requires three new 
parameters. The association constant for Atf3 protein with Atf3 
DNA {K\ i) was assumed to have an initial value of 0. 1 nM. The 
forward and reverse rates of the Atf3 protein for Atf3 DNA, Xi and 
X-i, were also required. X-\IX\ was assumed to have a similar 
value to All, but a slow rate of reversal giving a value of 
Xx = 1 xlO^ (Ms)^'. With the initial parameters, there was little 
effect on Atfi mRNA, Atf3 protein or Egri mRNA profiles 
(Figure 2, B— D, solid lines) compared with the original model 
(Figure IK). Decreasing the disassociation constant of Atf3 protein 
and Atf3 DNA by 10*" produced profiles that were more 
representative of the experimental data for Atf3 mRNA and 
protein (Figure 2, B and C, dashed lines), but acute inhibition of 
Egri mRNA expression was lost (Figure 2D, dashed line). 
Decreasing the ratio of (informed by a sensitivity analysis) 

2 orders of magnitude had no effect on the results. Increasing the 
ratio 2 orders of magnitude decreased the amount of Atf3 and 
increased the amount of Egri as expected (not shown). We 
conclude that Atf3 is unlikely to act in an autorepressive manner 
on its own promoter. If it did, the rate of accumtdation of Atf3 
mRNA/protein would be reduced such that it could then not 
significandy influence other gene promoters (in this case Egri) 
producing a redundant system. Thus, the Atf3 autorepressive loop 
is implausible. 

Model extension 2: Simultaneous phosphorylation of an 
alternative TF by RSKs inhibits Atf3 transcription 

Exposure of cardiomyocytes to BI-D1870, an inhibitor of RSKs 
[31], substantially enhances upregulation of Atf3 mRNA by ET-1 
(Figure 3A). Although BI-D1870 alone increase ERKl/2 phos- 
phorylation in cardiomyocytes, this over a similar time and to a 
lesser degree than that induced by ET- 1 and it does not affect the 
activation of ERKl/2 by ET-1 [9]. It is therefore unlikely that the 
enhancement of Atf3 mRNA expression induced by ET-1 seen in 
the presence of BI-D1870 is due to its elTects on phosphorylation of 
ERKl/2 and suggests that a negative feedback system down- 
stream of RSKs (or other kinases that may be inhibited by BI- 
D1870) moderates Atf3 mRNA expression. Interestingly, the 
effects of BI-D1870 over 2 h (Figure 3B) produced a profile for 
Atl3 mRNA expression that resembled the predicted profile for 
Atf3 mRNA in the modifications to the original model (Figure 1 K) 
suggesting that an input from RSKs might suffice to switch off Atf3 
transcription. We therefore extended the mathematical model to 
test the hypothesis that RSKs (as exemplar kinases) phosphorylate 



nuclear TFs (RTF) to negatively regulate Atf3 transcription 
(Figure 3C). Experimentally, maximal phosphorylation of RSKs 
occurs at 5-15 min, declining to ~25% the maximal level by 
~30 min (Figure 3D); this profile was modelled mathematically 
(Figure 3E, cyan line; see Text SI for details). The predicted 
profile for phospho-RTF binding to the Atf3 promoter is shown in 
Figure 3E (red line). Because different signals are appUed to 
different proteins, we assumed that the negative signal from 
phospho-RSK and phospho-RTF was competitive with the 
positive signal from phospho-ERKl/2 and phospho-TF. New 
parameters were required for the rate of RSK phosphorylation by 
phospho-ERKl/2 (Aig), the rate of RSK-P dephosphorylation ((//o), 
rates of association and disassociation of RTF for Atfi DNA {X^ 
and /t_2) and the rates of phosphorylation and dephosphorylation 
of RTF bound to Atf3 DNA (^lo and ^-lo). A fit-by-eye to the 
RSK data (Figure 3D) yielded A:9=lxl0^ (Ms)"' and 
(//o=:5.9xl0~* s~'. We assumed that = and /'._2 = A_i, 
k\(, = and k- lo = k^^, and that the total amount of RSKs is equal 
to the amount of ERKl/2 (i.e. = Eq). With the initial 
parameters, the profiles for Atf3 mRNA and protein remained 
similar to the original model (Figure IK), although the amount of 
each was reduced (Figure 3, F and G) and acute inhibition of Egri 
mRNA remained (Figure 3H). With the competition between 
phospho-TF and phospho-RTF for Atf3 DNA, increasing the rate 
of association of phospho-RTF with the Atf3 promoter may 
produce a profile that approaches that of the experimental data , 
but acute inhibition of Egri mRNA expression becomes lost (not 
shown). We conclude that the time delay between activation of 
ERKl/2 relative to RSKs is too small (1-2 min) to result in 
differential effects on mRNA expression over 0.5-4 h if both act 
directly on transcription from the same promoter. A greater time 
delay may be elicited by termination of transcription by an lEG 
encoding an inhibitory TF (ITF) and/or upregulation of a miRNA 
(or another factor) promoting Atf3 mRNA degradation and 
inhibiting Atf3 protein synthesis. 

Model extension 3: Upregulation of an inhibitory TF (ITF) 
that inhibits AtfB transcription 

We next tested the hypothesis that downregulation of Atf3 
mRNA expression (we retained the assumption that this is via 
RSKs) is mediated by upregulation of an ITF (Figure 4A), making 
the assumption that the negative effect of ITF was competitive 
with the positive signal from ERKs. Five further parameters are 
required for this model: the disassociation constant of ITF protein 
for Atf3 DNA (R'l^), the rate of ITF transcription (ky^, the rate of 
translation of ITF protein (^14), ITF mRNA degradation ((4), and 
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Figure 2. Model extension 1: Atf3 protein inhibits Atf3 transcription in an autorepressive loop. A, The conceptual model of Atf3 self 
regulation in which Atf3 protein suppresses the transcription of AtfB (dashed line). Symbols are as for Figure 1C. B, C and D, Atf3 mRNA (B, blue lines), 
AtfB protein (C, green lines) and Egri mRNA (D, pink lines) concentrations as predicted by the mathematical model using the initial parameter 
settings (solid lines) and with a 10^-fold increase in the association rate of AtfS for the Atf3 promoter (long dashed lines). In panels B and C, the 
graphs showing a 10^-fold increase in the association rate of Atf3 for the Atf3 promoter were also rescaled for comparison with initial parameter 
settings (short dashed lines). The predicted curves are superimposed on the experimental data (solid black lines, right axes) to which a spline curve 
was fitted. 

doi:1 0.1 371 /journal.pcbi.1 003597.g002 



the rate of ITF protein degradation (rf-,). On the basis that ITF is 
an lEG, we assumed that the disassociation constant of ITF 
protein was equal to that of TF for Atf3 DNA (A'ls = Kiq), rates of 
transcription and translation of ITF are similar to those for AtfS 



(i.e. kj:j = k(i and k^4 = kji), and the rates of degradation of each 
mRNA and protein are similar (i.e. = rfj and (4 = d^) producing 
the profiles for ITF mRNA and protein shown in Figure 4B. With 
initial conditions in which ITF binding to AtfS DNA is similar to 
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Figure 3. Model extension 2: Simultaneous phosphorylation of an alternative transcription factor by RSKs inhibits Atf3 
transcription. A and B, Experimental data for the effects of inhibiting MKK1/2 with PDl 84352 (2 \xM) or inhibiting RSKs with BI-D1870 (10 jiM) on 
AtfB mRNA expression in unstimulated cells or in cardiomyocytes exposed to ET-1 (100 nM, 1 h; A) or for the times shown (B). Results are expressed 
relative to controls as means ± SEM for n = 3 (A) or 4 (B) independent myocyte preparations. C, The conceptual model was developed to include 
phosphorylation of RSKs (RSKP) that phosphorylate a RSK-sensitive TF (RTF) to produce phosphorylated RTF (RTFP) that binds competitively to the 
AtfB promoter to suppress AtfB transcription. Symbols are as for Figure 1 C. D, Experimental data for RSK phosphorylation in cardiomyocytes exposed 
to ET-1 for the times shown. Results are means for n = 3 independent experiments and are adapted from the data in [9]. E, Predicted profiles from the 
model for RSK phosphorylation (cyan) and RTFP binding to the AtfB promoter. F, G and H, AtfB mRNA (F, blue line), AtfB protein (G, green line) and 
Egri mRNA (H, pink line) concentrations as predicted by the mathematical model using the initial parameter settings The predicted curves are 
superimposed on the experimental data (solid black lines, right axes) to which a spline curve was fitted. 
doi:1 0.1 B71 /journal.pcbi.l 00B597.g00B 



Atfi binding to tlie Egrl promoter (Figure 4, C and D, solid lines), 
there was litde change in the profiles compared with the original 
model (Figure IK). Increasing the rate of binding of ITF to Atf3 
DNA by 10^ produced profiles that were more similar to the 
experimental data (Figure 4, C and D, dashed lines), but inhibition 
of Egrl transcription was reduced (Figure 4E, dashed line). 

Model extension 4: Upregulation of one or more miRNAs 
that regulate AtfB mRNA and/or protein expression 

Transcriptomics studies using AfFymetrix exon arrays [9] 
identified several miRNAs that are regulated in cardiomyocytes 
by ET-1 within 1 h and whose upregulation is inhibited by BI- 
D1870 (Table 1; since we have not performed miRNA profiling of 
cardiomyocytes, this is unlikely to be a full representation of 
miRNAs in this category). Atf3 expression may be regulated by 
miRNAs (e.g. Mir663 suppresses Atf3 mRNA expression in 
endothelial cells subjected to shear stress [32]). We therefore tested 
the hypothesis that upregulation of one or more miRNAs (we 
retained the assumption that this is downstream of RSKs as 
exemplar kinases) is required to downregulate AtfS mRNA 
expression (Figure 5A). This extension to the model introduces 
eight new parameters. With little/no published information on 
rates of miRNA turnover, we informed the model as follows. The 
rate of phosphorylation/dephosphorylation for RTF binding to 
miDNA was assumed to be equal to the rate of phosphorylation/ 
dephosphorylation of TFs bound to Egrl and Atf3 DNA by ERK- 
P (kif, = k-i; k-i(, = k-'ij. The rate of miDNA transcription was 
assumed to be the same as that of Atf3 transcription [kn = k^) and 
the rate at which premiRNA is converted to miRNA (kig) was 
assumed to be double the rate of Atf3 translation (0.5/s). The rate 
at which miRNA associates with AtfB mRNA {kig) was assumed to 
be equivalent to the rate set for AtfB protein binding to Egrl DNA 
(1 xlO' (Ms) ') and the reverse rate (^-19) was set at 5x10 ^/s, 
the value that gave the best qualitative fit to the experimental data. 
The rates of decay of premiRNA (df,) and miRNA (rfy) were 
assumed to be ~24 h in line with recently published data [33,34]. 
The rate of decay of the complex miRNA- mRNAAtfB (dg) was 
initially assumed to be the same as that of AtfB mRNA {dg = d^). 
Mathematical modelling for the production of miRNA is shown in 
Figure 5B. 

miRNAs may increase mRNA degradation (within a complex 
with Ago2) or silence mRNA expression (within Agol, AgoB or 
Ago4). Computational modelling has been used to predict/ 
confirm that the increase in mRNA degradation induced by 
miRNAs lies within the range 1.3- to 6.4-fold P)5]. We first 
assumed that potential miRNA(s) had a negligible effect on the 
rate of degradation (i.e. the rate of degradation of the miRNA/ 
mRNA complex was similar to the free mRNA). In this case, the 
total AtfB mRNA (i.e. free AtfB mRNA + miRNA/mRNA 
complex) dynamical profile did not match the experimental data 
(Figure 5G, solid line). However, since the miRNA/mRNA 
complex is assumed to be translationally incompetent, the protein 
profile was significantiy affected and approximated to the 



experimental profile (Figure 5D, solid line). If we assumed that 
potential miRNA(s) increased the rate of degradation of Atf3 
mRNA from die miRNA/mRNA complex 2-fold, the total Atf3 
mRNA profile approached that of the experimental data 
(Figure 5C, long dashed line). Increasing the rate of degradation 
4-fold gave a profile that approximated to the experimental data 
(Figure 5C, short dashed line). The protein profile was relatively 
unaffected by increasing the rate of degradation of the miRNA/ 
mRNA complex (Figure 5D), indicating that miRNAs have a 
dominant effect on protein expression irrespective of their effect on 
mRNA expression. In all cases, the predicted Egr 1 mRNA profile 
approximated to the experimental data over the 2 h period studied 
(Figure 5E), although subsequendy Egrl mRNA levels were 
predicted to increase as AtfB protein is lost. 

Combining model extensions 3 and 4: Upregulation of 
both an ITF and miRNA(s) is required for inhibition of 
AtfB mRNA and protein 

We combined the two systems for ITF and miRNA inhibition of 
AtfB to investigate whether an ITF may influence the profiles in 
the presence of miRNA(s) (Figure 6A). Profiles for total Atf3 
mRNA (Figure 6B), AtfB protein (Figure 6C) and Egrl mRNA 
(Figure 6D) were similar to those for model extension 4 (Figure 5, 
C~E) indicating that miRNAs alone could suffice to switch off AtfB 
mRNA and protein expression. 

Discussion 

In considering mechanisms of regulation of lEG expression, the 
emphasis is usually on the signals and transcriptional regulators 
required for upregulation. Here, we explored possible mechanisms 
associated with downregulation of lEG expression and termination 
of the response. We focused on AtfB as an lEG in which there is 
great interest because of its association with inflammation, cancer 
and cardiac dysfunction [11-14]. 

The report that Atfi can bind to its own promoter to inhibit 
transcription [22] is, at first sight, attractive. The experimental 
approaches used remain the most appropriate, but all such studies 
are potentially subject to artefacts resulting from overexpression or 
elimination of a single transcription factor in isolation. Developing 
a mathematical model provides an alternative means of testing the 
hypothesis. In this case, the model demonstrated that such a 
negative feedback loop is implausible in the context of our 
cardiomyocyte system (Figure 2). This raises the question of 
whether any such self-regulating transcriptional system can 
operate. For an inhibitory TF, the model prediction is that a 
redundant feedback loop is most likely to be generated in which 
the TF serves to regulate only its own expression, restricting the 
possibility that it accumulates to a sufficient degree to operate in 
an efficient inhibitory loop with another gene. The scenario would 
differ in the case of a positive TF driving expression of second 
phase genes. Here, self-regulation can provide a means of 
moderating downstream gene expression. 
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Figure 4. Model extension 3: Upregulation of an ITF that inhibits Atf3 transcription. (A) The conceptual model of the ITF inhibitory 
pathway for Atf3 transcriptional regulation (dashed lines). Symbols are as for Figure 1C. B, ITF mRNA (left axis, solid line) and ITF protein (right axis, 
dashed line) concentrations as predicted by the mathematical model. C, D and E, Atf3 mRNA (C, blue lines), Atf3 protein (D, green lines) and Egri 
mRNA (E, pink lines) concentrations as predicted by the mathematical model using the initial parameter settings (solid lines) and with a 10^-fold 
increase in the association rate of ITF for the Atf3 promoter (long dashed lines). In panels C and D, the graphs showing a 10^-fold increase in the 
association rate of ITF for the Atf3 promoter were also rescaled for comparison with initial parameter settings (short dashed lines). The predicted 
curves are superimposed on the experimental data (solid black lines, right axes) to which a spline curve was fitted. 
doi:10.1371/journal.pcbi.1003597.g004 

Our experimental data suggest that, whereas ERKl/2 them- important in terminating Atf3 mRNA and, presumably, protein 
selves promote upregulation of Atf3, downstream activation expression (Figure 3, A and B). This resulted in the development of 
of RSKs (or other BI-D1870-sensitive kinases) is potentially three further extensions to the model. Whilst we have included 
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biologically appropriate mechanisms at each model revision one 
outcome of this has been an increase in the signalling cascade 
complexity. This in itself has led to a more robust signalling 
response. 

The initial concept, that delayed activation of RSKs may result 
in delayed activation of a pre-existing inhibitory TF for signal 
termination, appeared attractive. However, model extension 2 
demonstrated that direct activation of the putative RSK-respon- 
sive RTF was almost simultaneous with activation of the ERKl /2- 
responsive TF and simply titrated down the levels of Atf3 mRNA 
rather than altering the profile. Thus, the delay between activation 
of ERKl /2 and activation of RSKs is insufficient for this scenario 
to be viable and a greater delay is required before initiation of the 
negative signal. Interestingly, since a delay is necessary, this also 
negates the possibihty of other acute events (e.g. histone 
deacetylation) being responsible for terminating Atf3 mRNA 
expression. The timing strongly suggests that a second transcrip- 
tional event is required, either to produce an alternative inhibitory 
TF (termed ITF, model extension 3) and/ or one or more miRNAs 
(model extension 4). Here, developing first model extensions 3 and 
4 independently and then combining the two strongly suggested 
that an ITF would not affect Atf3 mRNA and protein expression 
whereas miRNA(s) most probably would suffice to generate the 
profiles seen experimentally. Interestingly, the miRNA effect was 
particularly dominant on the Atf3 protein profile. The model also 
demonstrated that a modest increase in the rate of degradation of 
Atf3 mRNA by the miRNA within the range predicted by other 
computational models [35,36] could suffice for the Atf3 mRNA 
profiles to approximate to the experimental data. Further studies 
of cardiomyocyte miRNAs are clearly required. 

Materials and Methods 

Quantitative PCR and microarray data 

Neonatal rat cardiomyocytes were prepared as previously 
described [37] and exposed to 2 |xM PD 184352 (to inhibit the 
ERKl/2 cascade [38]) or 10 nM BI-D1870 (to inhibit RSKs [31]) 
for 70 min, to 100 nM ET-1 for 1 h or to inhibitor for 10 min 
prior to addition of ET-1 (1 h). PD184352 and BI-D1870 were 
from Enzo Life Sciences and were dissolved in DMSO. ET-1 was 
from Bachem UK. RNA was extracted and quantitative PCR 
performed as described in [10]. Affymetrrx microarray data for 
cardiomyocytes exposed to ET-1 in the absence or presence of 
PD184352 or BI-D1870 were published in [9]. The data are 
available from ArrayExpress (accession nos. E-MIMR-3, E- 
MIMR-37, E-MEXP-3393, E-MEXP-3394, E-MEXP-3678, E- 
MEXP-3679). Probeset sequences that were significantiy upregu- 
lated and not assigned to a protein-coding gene were analysed by 
BLAST search against the rat and mouse genomes (www.ncbi. 
nlm.nih.gov/ genome/ seq/BlastGen/BlastGen.cgi?taxid = 10116; 
www.ncbi.nlm.nih.gov/ genome/ seq/BlastGen/BlastGen.cgi?taxid 
= 10090; cross-species megaBLAST) to identify those which 
recognise miRNAs. 

Data analysis and curve-fitting to experimental data used 
GraphPad Prism. 

Mathematical modelling 

FuU details of the mathematical modelling are provided in Text 
SI. 

Sensitivity analysis. A thorough local sensitivity analysis of 
each of the mathematical models was undertaken in which each 
parameter value was varied 10- and 100-fold above and below its 
originally determined value in fitting to the experimental data 
(phosphorylated ERK, phosphorylated RSK, phosphorylated 
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Figure 5. Model extension 4: Upreguiation of one or more miRNAs inhibits Atf3 mRNA and/or protein. A, The conceptual model of the 
mlRNA inhibitory pathway for Atf3 transcriptional regulation. Symbols are as for Figure 1C. B, The predicted dynamical variation in mlRNA 
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lines) and Egri mRNA (E, pink lines) when the rate of degradation of AtfS mRNA within the miRNA.mRNA complex is the same as for degradation of 
free AtfS mRNA (solid lines), is increased 2-fold (long dashed lines) or is increased 4-fold (short dashed lines). The predicted curves are superimposed 
on the experimental data (solid black lines, right axes) to which a spline curve was fitted. 
doi:1 0.1 S71 /journal.pcbi.1 00S597.g005 
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CREB, phosphorylated c-Jun, Atf3 mRNA, Atf3 and Egrl 
mRNA). Varying the parameter values associated with the 
extension of our original mathematical model whereby phosphor- 
ylation of transcription factors for Atf3 and Egr 1 were explicidy 
included, did not give any unexpected results, for instance 
variation in the rate of MKK 1/2 or ERK phosphorylation did 
not alfect the fit to the data. The fit to Egrl mRNA data was found 
to be most sensitive to those parameters which decrease the 
amount of available Atf3 (decreasing Atf3 transcription, decreasing 
AtfS translation, increasing Atf3 mRNA and Atf3 degradation) 
and the rate at which AtfiS supresses Egrl mRNA. 

A sensitivity analysis of the two models in which either Atf3 
protein or RSKs inhibit Atf3 transcription showed that no increase 
or decrease in parameters associated with each model improved 
the fits shown in Figures 2 and 3, respectively. The fit of each 
model was particularly sensitive to how much Atf3 was produced 
(for instance increasing the rate of RTF association for Atf3 DNA 
l()-f()ld simply decreases the amount of AtfS such that the amount 
of Egrl mRNA increases, removing the transiently observed 
decrease around 2 h). 

We found that the fits between the remaining three models (ITF 
inhibition of Atf3, miRNA regulation of Atf3 mRNA and a 
combination of the two) and the experimental data were 
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particularly robust to variations in parameters associated with 
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reverse phosphorylation of the ITF bound transcription factor 
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of the model solutions; any increase or decrease in ITF simply 
affected the Atf3 concentration. Similar results were found for the 
miRNA model extension and the combined extension of ITF and 
miRNA upregulation. 
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